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Abstract 

The  paper  is  concerned  with  the  decay  of  the  energy  of  disturbances 
which  are  propagated  according  to  the  wave  equation  with  variable  index 
of  refraction  in  the  exterior  of  a  finite  star-shaped  reflecting  body.   It 
is  shown  that  the  energy  of  the  disturbance  decays  like  some  power  of  t 
Certain  conditions  of  growth  and  continuity  are  made  on  the  index  in  order 
to  insiire  some  decay  factor. 

The  energy  decay  is  obtained  by  estimating  the  solution  of  an  integral 
equation  which  results  when  one  applies  the  Friedrichs  '  "A-B-C  method"  to 
the  modified  wave  equation  operator.   Using  the  energy  estimate  with  other 
familiar  estimates^  one  obtains  a  rate  of  decay  for  the  disturbance  itself. 


^v 


Introduction 

Let  us  consider  disturbances  which  are  propagated  In  three-dimensional 
space  according  to  the  modified  wave  operator 

u   +  u   +  u   -  c  '^fx.y^z)  u,  ,  =  0. 
XX    yy    zz      ^  •'''  ^  '   tt 

If ^  initially,  the  disturbance  were  confined  in  some  finite  region  in  the 
exterior  of  a  finite  smooth  reflecting  body,  one  would  expect  that  after 
a  certain  fixed  time  most  of  the  distiorbance  would  have  propagated  to  Infinity. 
In  [1]  it  is  shown  with  certain  restrictions  on  the  body  that  this  is  the 
case  for  c  =  1,  and  the  disturbance  at  a  point  dies  out  exponentially. 

For  the  non-homogeneous  wave  equation  one  also  finds  that  if 
u  satisfies  Initial  data  of  compact  support  and  the  equation 

/     ^  i(jJt 
u   +  u   +  u   -  u^_^  =  g(x,v,z)e 
XX    yy    zz    tt   ^^  '"  '  ' 

in  the  exterior  of  a  smooth,  finite,  star- shaped  body  on  which 

u(x,y,z,t)  =0  t  >  0, 

then  ue     tends  exponentially  in  time  to  v(x,y,z)  at  every  fixed  point 
exterior  to  the  body.   Here  v  satisfies 

2 
V   +v   +v   +wv=s 
XX    yy    zz 

exterior  to  the  body, 

V  =  0 

on  the  body  and  v  satisfies  the  Sommerfeld  radiation  condition  at  Infinitv. 


This  result  is  known  as  the  limiting  amplitude  principle.   The  proof  is 
based  on  Kirchoff fe  theorem  and  on  the  previous  energy  estimate. 

One  expects  that  both  the  energy  decay  and  the  limiting  amplitude 
principle  are  also  true  with  a  variable  speed  of  propagation.   Here  an 
estimate  will  be  found  for  the  energy  decay  in  any  fixed  region.   A  proof 
of  the  limiting  amplitude  principle  similar  in  nature  to  [2]  is  not  possible 
since  Kirchoff 's  representation  theorem  is  no  longer  valid  in  an  inhomogene- 
ous  medium.   But^  it  will  be  shown  that  all  solutions  to  the  initial 
boundary  value  problem  for  the  homogeneous  wave  equation  with  variable 
index  of  refraction  decay  in  time  at  any  fixed  point  exterior  to  the  body. 

One  may  expect  some  restrictions  on  the  index  of  refraction  for, 
if  the  index  oscillated  very  rapidly,  the  energy  might  be  trapped  between 
two  minima  of  the  speed  of  propagation.   However,  in  the  case  of  the  earth 
and  its  atmosphere,  a  good  approximation  for  the  behavior  of  the  index  of 
refraction  is : 

c(r)  =  c   -  c,  e 
^  '    o    1 

where  c  >  c^  >  0  and  p  >  0. 
o    1 

This  case  will  be  covered  by  the  hypotheses  on  the  index,  and  the  energy 

does  go  to  infinity. 

The  estimates  are  obtained  by  the  methods  used  in  [2] j  i.e.,  the 

energy  estimate  will  be  obtained  by  a  variation  of  the  A-  B  -  C  method 

of  Friedrichs.   This  consists  in  multiplying  the  differential  operator 

_2 
A  u  -  c   u,  ,  =  0 


by  the  operator 


Au  +  B'Vu  +  Cu 


-  3  - 


and  integrating  over  the  space  exterior  to  the  body.   By  strategically 
choosing  the  A^  B,   and  C,  and  by  applying  the  divergence  theorem  to  the 
resulting  integral^  one  obtains  an  integral  equation.   Upon  estimation  of 
certain  terms  in  the  integral  equation,  one  finds  that  the  energy  in  any 
fixed  finite  region  exterior  to  the  body  decays  in  time  like  some  power 
of  t"'^. 


The  main  theorem  to  be  proved  is: 

Theorem  1:   'Let  u  be  a  smooth  solution  of 

(l)         u   +  u   +  u   -  c  {x,y,z)u..    =  0 
^   '  XX    yy    zz      ''  '-^  ■'  '  tt 

in  the  exterior  of  a  smooth,  star-shaped,  finite  three-dimensional  body  B 
on  which  the  boimdary  condition 


(2)  u  =  0 

is  satisfied.   Let  u  satisfy  the  initial  conditions 

(3)  u(x,y,z,o)  =  f(x,y,z) 

(^)         u^(x,y,z,o)  =  g(x,y,z) 

where  f  and  g  are  both  twice  continuously  differentiable  and  of  compact 
support.   Suppose  c(x,y,z)  satisfies  the  following  conditions: 

(5)     I     c(x,y,z)  is  continuously  differentiable 

II     c(x,y,z)  is  a  function  of  the  distance  from  the  origin  of  the 
coordinate  system 
III     min  c(x,y,z)   =  7  >  o 
IV     |c|  <  min  (^   ,     -\ 


for  all  r 


.r~2         2    2  ^  ^ 
/x  +y  +z  >  0. 


Then 


(6)   t^"  ^/   (|V  (ru)  f  +   c"^(ru^)^)  sin  9  drdcpdQ  <  M 


R 


-  5  - 


where  M  is  some  constant  depending  upon  c,  f^  g^  R  and 


(7)     a  =  max    :^  I  |£  I  <  5  <  2  . 


Before  proceeding  with  the  proofs  a  few  remarks  will  be  made  on  the 
conditions  on  the  index  of  refraction.   Condition  I  is  needed  for  the  appli- 
cation of  the  divergence  theorem  to  the  differential  operator.   Condition  II 
was  assiimed  because  of  the  simplification  involved  in  the  integration  and 
also  because  of  its  good  approximation  to  physical  considerations.   It  seems 
possible  that  angular  dependence  could  be  introduced,  provided  certain 
boundedness  conditions  are  made  on  the  index  and  its  derivatives.   Condition 
III  is  needed  to  maintain  the  hyperbolic  character  of  the  differential  equation. 
Condition  IV  puts  the  greatest  restriction  on  the  index.   In  effect,  it  says 
that  for  large  r  the  index  should  approach  a  constant  like  0(r  ) .   That  this 
condition  is  reasonable  can  be  seen  from  the  paper  by  Miranker  [ 3]  in  which 
he  shows  that  the  Rellich  uniqueness  theorem  is  valid  for  the  reduced  wave 
equation  with  variable  index  of  refraction,  provided  the  index  approaches  a 
constant  for  large  r  like  0(r  ) .   Thus  this  condition  is  a  natural  one  if 
a  limiting  amplitude  principle  is  to  hold  for  the  variable  case.   That  |c| 
should  be  less  than  7/2r  might  be  improved;  however,  it  was  the  simplest 
necessary  condition  in  order  to  get  a  decay  in  the  energy.   If,  for  example, 
we  consider  the  earth's  atmosphere, 

c(r)  =  1  -  (1  -  7)e-'^'', 
we  obtain  the  following  restriction  on  7.   [We  assiitne  that  c(r)  has  been 


-  6 


normalized  to  be  one  at  infinity.]   By  (5) 
|c(r)|  =  (1  -  7)|3e"'^''<  7/2r 


or 


2(1-7)  P7"^  <  eP^'r"^. 


The  minimum  on  the  right  side  is  obtained  for  r  >  0  when  r  =  3   .  We  obtain 


2(1  -  7)  7""""  <  e 


or 


2(2  +  e)"^  <  7. 
Thus  the  index  may  change  by  no  more  than  e(2  +  e)    of  its  maximum  value 
in  order  for  our  result  to  hold.  We  now  proceed  with  the  proof  of  the 
theorem. 

Proof:   Choose  a  spherical  coordinate  system  so  that  the  body  B  is  star- 
shaped  with  respect  to  the  origin,  i.e.,  r  <  0,  where  r  is  the  radial 
component  of  the  ;init  normal  pointing  into  the  body.  Let  D,  be  the  domain 
of  influence  of  the  initial  data  at  time  t.  Writing  equation  (l)  in 
spherical  coordinates,  one  obtains 


(Q)      (r%sin  9)^       (n^   sin  9)^      ^cpcp        \t    _  _ 

? "*"     2  "*"  2   2   "   2 

r  sin  9  r  sin  9       r  sin  9    c  (r) 


Multiplying  (8)  by  the  operator 

(9)     Au  +  Bu  +  Cu_^ 
^^'  r     t 

and  integrating  over  the  entire  domain  of  influence,  one  obtains 
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r   1   r 


/(r  u  sinO)   (V^'^^^g   ^cpcp 
(10)  J         J         (AU  .  BU^.  CU^)(  -.^i— ^.  -.^-^  .   -« 

o     D 


tt 


2        '22 
r'"sin0      r  sinO      r  sin  9 


t 


Now  set  ru  =  w  and  choose 

(11)  A  =  2t 

(12)  B  =  2rt 


(13) 


C  =  t  + 


r  c^(0 
o 


d^ 


dvdt=  0, 


For  convenience  the  explicit  form  of  B  and  C  will  not  be  exhibited;  however, 
we  will  exploit  their  form  in  the  following  calculations.   After  the  change 
of  variables,  (lO)  becomes 

(1^)    ^ 

(w^sin  e).   w^ 


'     /   (Bw  +Cw  )  (sinOw  + 


9    cpcp       tt     \   dvdt  ^ 


o     D, 


2     ■   2  .  ^     2,  .    /  2 
r        r  sin9    c  [r)        I   r 


Using  the  fact  that 

2  2 

2(Bw  w   )  =  (Bw  )    -  B  w 

2       2 
and  2  (Cw^w  )  =  2(Cw  w_^)   -  2  C  w  w ,  -  (Cw  ),+  n  w   ,   etc., 
^   t  rr        r  t  r      r  r  t      r  t   t  r   ' 

one  obtains  for  (lU) 


o         D, 


Bsin9w^  +  2Csin9w  w_^-  -^"« 
r  r  t        2 


Bsin9  2 


Bsin9w 


9 


r      sm   9     c    (rj  r 


/  2  2  2  \ 

pon-r^chr^    ,    Cw  „     ,    2sin9Bw,w      ,    Csin9w^        Csin9w_    \ 
Csin9w^  +__cp    +  _  t  r  +     ^-  t  +  -^  9 

\  r  sin9  c    (r)  c    (rj  r  /   t 

■2Bw  w  2Cw  w^  \        /2sin9Bw^w  2Csin9w^w^ 


2  ■      2 

r  sin  9       r  sin  9 


9 


1    2^c(r) 
+     Uw  rt   4   ' 

t     c^(r) 


sin9 


drd9d9dt  =   0. 
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Applying  the  divergence  theorem  to  this  equation  one  obtains 


(15) 


2Cw  w^  +  B(2w  - 
r  t       r 


IVwl^  + 


w, 


c^(r) 


[ 


2^  w^w^  .C(|Vv|2.!|^   ) 
c  (r)       .  c  (r)   J 


2   2        2   2 
'vr  sin  0      r  sin  9   /  9 


V 


/  2B  WW     2C  w^w^  \  ^ 
+  ( —  9  r  +  -^  9  t  j  9^ 


sin9dSdt 


>t,  n 


/    y-^  c(r)  w^sin9draW9dt  =  0, 


D, 


c^(r) 


where  r  ^  t  ^  cp  .  and  9  are  the  various  components  of  the  unit  normal 

VVV  V 

pointing  out  of  the  four-dimensional  time-space  region.   The  surface 

integral  in  (I5)  may  be  subdivided  into  three  integrals  :  one  over  the  body 

B  from  time  t  =  0  to  t  =  t-,  j  the  second  over  the  initial  data  up  to  the 

radius  r  =  a  ;   the  third  over  the  region  D,  .   There  is  no  contribution 
o  t-, 

for  large  values  of  r  since  the  solution  is  identically  zero  outside  the 
domain  of  influence  of  the  initial  data  which  are  of  compact  support. 
The  integral  over  the  body,  I  ,   is  given  by 


(16) 


^1^ 


o    B 

it 


B(2w2  -  |Vw|2)  r  .  ^  Vr^v   .  2  5  Vr^^ 
r    '   '    V   2  .  2^  2 

r  sm  9         r 


sin  QdSdt 


1  /  2t   /dw\       -,„^^ 
r   \dn/     V 

B 


-  9  - 


where  we  have  made  use  of  (2) ,    (l2)  and  the  fact  that  t  =  0  on  B.   Since 
B  is  star-shaped  r   <  0  and  therefore 

The  integral  over  the  initial  data  is 


(18)    I^  =   /  (  /   -^    a  A       I  |V  wl^  +  ^    1  sin9drdcpdQ 


'       i^o    ='(5)    'IV'     Ar) 


We  make  use  of  (l3)  to  obtain  (l8).   By  (5)^  (^) ^  and  (5)  I   is  bounded. 

The  integral  over  the  region  D  ,   1    ,   is  written  as 

*1        ^ 


(19)  I5 


=   "/  r^^'t^r  +  c(|V  w|2  +  ^\   1  sin0drdq,d0. 
^     Lc    (r)  \^  c    (r)y    -" 

*1 


Combining  (16) ,    (18) ,   and  (19)  and  using  (l7)  ,   one  obtains 

\ 

-  -^  Urt  ^4^     wj   sin  9  drdcpdOdt        <  I  ^  * 

t^ 

Since  j        /    ^-   |c(r)|  |Vw|   sin  ©drdcpd9dt   is  positive  we  can  subtract 


o   \       =-(^) 


it  from  the  left  side  of  (20)  without  destroying  the  inequality. 


10  - 


Consequently  we  have 

r 


2./'^  i^  a.    )    f  lv.|2  .  "^  \.^^t. 


^P  ^^     c  (o        V  c  (^)/  c  (^) 


*1  tdt    r    ^^i44^         [  |Vw|2  +     !L      1      sin    Qdrdcpde  <  I 


cirj     \  '       •  „2/  \ I  2 


Considering  the  first  integral  in  (21)  ,   we  integrate  /    f-"   <i).   by 
parts  to  obtain  for  this  term  o 

(22)  f  l(t^  +  £-       )       Avw|2  +  ^y^  WW  V  sinOdrdcpdQ 

J^  [V       c2(r)/   ^^       c2(r)y   c2(r)   ^  ^ 


*1 


+    /  1     /  2^ 

D,  \  ^o 

^1 


^2       \ 
■2     £lLl_    d     1  IVwi^  +  -^ —    1       sin9drdcpd0 


\)     V  '       c=^M, 


Applying  Schwarz  's  inequality  to  the  term  i+rt  w  w  c   (r)  in  (22)  , 

X  Xf  j/ 

one  has 

(25)     |lfrt^w^w^c'^(r)|  <  2rt^c"^(r)  (w2  +  c'^w^)  • 

Subtracting  the  positive  term 

D^  '    \  r  sm  0      r     / 

from  the  left  side  of  (21) ,  and  applying  (23)  to  the  first  integrand  does 
not  alter  the  inequality,  and  one  obtains 


11 


r       \    2       /,„    ,2         -t 


{2k)        I        {\-^  U"""'      '^)1        -^^^^^^'^ 


^H 


^1  tdt    r    ^^1^1  [  |Vw|2  +  4;i_  ]      sin  9  drdcpde 

1^ 


<  I   _  /     [  /     "S^"\'a/  d.  I   (  |Vw|^  +  -^     \    sin9drdcpd9 


Wow  if  we  suppose  that  c(r)  is  a  monotonia  increasing  function,  the 
second  term  on  the  right  side  of  {2k)   would  be  strictly  positive  and  we  could 
therefore  neglect  this  term  since  the  left  side  would  then  surely  "be  less 
than  I  .   However;,  if  c(r)  is  not  monotonia  increasing  we  proceed  as  follows. 
In  (10)  set  A  =  B  =  0  and  C  =  1.   One  thus  has 


w  isinOw  ^+  —^ 9  +^li_--^— -     -2-dt  =  0 

\  r         r  sm  9    c  (r)    /   r 

o    D^     V  ^     / 


Integrating  as  before  one   obtains 


^\    r      {    ,  ,         .      ^       ^   (ytV  ,      (vQsin9w^ 

I  r   sin  9  r 


r  sin9         c    (r)  r  /  t  1    r 
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Applying  the  divergence  theorem,  there  resiilts  the  following  surface  integral 


/         2w^^    2w^v^9^ 
\  r  sm  0    r 


sin  0 


IVwl^  +  -i —  It  sin  0 
'       c^Cr)/   - 


dSdt  =  0. 


As  before, the  surface  integral  is  taken  over  the  body  B,  the  initial  domain 
and  the  domain  of  influence  at  t  =  t^ .   Since  w  =  0  on  B,  w  =  0  on  B. 
Therefore  the  integral  reduces  to  the  following  identity 


(25)     /     (l 


,2         ""t 
Vw       +  —T —     )  sin  0  drd9d0 

c^(r) 


D 


IVwl^  +  -^     )sin  0  drdcpd0. 
c2(r)  / 


This  identity  expresses  the  notion  of  consei^ation  of  energy. 
Wow  consider  the  last  term  in  {2k). 


(26) 


D 


t   \  o 


|Vw|2  +  ""^ 


c^(r) 


sin  0drdcpd0 


< 


max 


r  <3:<oo 

<-    o  "^  r 


J-  c5(0        - 


|Vw|^  +  -^     pin  0  drdq)d0 
li'    D    \  c^(r)/ 
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The  last  integral  is  finite  by  conditions  (j)  and  (k) .      By  condition  (5) 
IV,  one  has 


Max 

■  <r<co 

o 

P""          2 
o 

f2) 

r\     +  ^  /   1/5  _L 


< 


1      i.^    .  2^  ^2l^ 


Therefore  the  left  side  of  (26)  is  bounded  by  some  finite  number  which  will 
be  called  T    . 

Choose  a  sphere  0  which  is  contained  in  D  .   Let  0  be  centered  at 
the  origin  and  have  radius  6  yt ,   where  0  <  ^  <  <  1.   Using  (7)  and  (26)  in 
{2h)      one  obtains 


2         ""t      \ 

+     -p-       sin  0  drdcpdO 

*1 
t.  /  2 

P    1  r      I  o  ^f 

V  w  I      +  — p-  I     sinQ  drdcpde 

°t    ^ 


(27)      t^   (1  -   s   f  J      hVw| 
t-, 


p 

<  I^  +  r     +   /     ^     tdt    r  ^^1^1  [    |Vw|2  +     -|-J  sin  9  drdcpdO 

Xi  u 


Using  (5)  IV  and  (25)  on  the  last  term  in  (27)  one  obtains 


11^ 


tdt 


Uriel        /,„    ,2       ""t 


sin  9  drdrpdQ 


< 


L  ^D 


|Vw|^  +  -^   I    sin  0  drdcpdO 


2  tdt  +    /  -^^ 


< 


,2       ^t 


L^D 


Vw  I      +  — ^  I  sin  0  drdcpd9 


"2      -3-14.-1 
T^  +  2e  7  T^ 


o 


-1-1   -1  ^/5 
where  T,  =  e  7   (27   )    .   The  last  expression  is  finite  "by  (5)  and  {k) 

and  is  bounded  by  some  constant  T    . 


Now  if  we  set 
2 


/   (I'^^i^  +  ^  )  sin  Q  drdq)d9  =  E  (O,), 


one  obtains  the  following  inequality  from  (27) 


(28)  t^  (1  -  e  f   E(0^  ) 


a   /    t  E  (O^)dt  <  I2  ^  -^1  ^  ^2  "  -^J" 


This  can  be  rewritten  as 

(29)   t^  E  (0^  )  -  a^  /    t  E  (0^)  dt  < 


where  cL,=a(l-6) 


^=  ^3  (1  -  e  )■ 


Set  t  E(0, )  =  G(t)  in  (29)  which  becomes  now 


-  15 


t  G(t)  -  a^  /  G  (T) 


dT  <  r,    or 

4- 


(50)     t    ^   |^(t    1  j     G  (T)dTJ  <r^. 

Integrating  (jO)  from  t  =  1  to  t  =  t^  one  obtains 

-QL,  o   t        r,    r,      ^1 

(31)   t^^J  G  (T)dT  <  -1-   -  -it^ai    .  j    T  E(0^)dT  . 

^  o  T.     1       "^  o 

[Note:  The  choice  of  the  initial  point  is  quite  arbitrary  and  t  =  1  was 
chosen  just  for  convenience.] 

But  by  (25)  the  last  term  in  (jl)  is  finite  and  (51)  becomes 


(32)    /'^^    G  (T)dT<  (^   +r^]  t""^ 


where  P^ 

5 


:/     T  E   (O^)d'r. 


Ret\irning  to   (29)  one  obtains 

t?  E(0^  )  <  r,   +  a,  I — -    +  r^  I  t  ■"" 

1  ^  tj_'         h         ll    (x^           5y 

or 


For  large  t    ,  this  becomes 


E(0^    )   <     M  t"^^""^!^ 


^1 

E(0,    )   <  Mt"^^"^^ 
^1 

0  <  e  «  1. 
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By  letting  t  become  large  enoiigh^  any  finite  fixed  region  wiUL 

eventually  be  included  in  0  .   Thus  for  large  t  ,  R  c^  0   and  by  the 

^1  ^  *1 

definition  of  E  (O  ),  one  obtains  the  result  as  required  in  (6).  Thus  if 
^1 

max      kr    \ . . 

r  <r<po    c  '  ' 

o 

is  less  than  2,   one  can  obtain  a  rate  of  decay  for  the  energy  in  any 
finite  region. 


IT  - 


Appendix 

It  is  now  possible  to  say  something  about  the  decay  of  the  solution 
u  of  (l)  -  (5) .  We  have  the  following: 

I 
Theorem  2:  If  u  is  a  smooth  solution  of  (l)  -  (5)  in  the  exterior  of  a  smooth, 

star-shaped  finite  body  B,  then 


(51^)    |u|  <M2  t-(^-^/2) 


The  proof  makes  use  of  two  basic  lemmas  which  will  only  be  stated  as 
they  are  quite  familiar  and  appear  in  [2] . 

Lemma  1.   Any  solution  u  of  (l)  and  (5)  satisfies  the  inequality 


(35)   |u(X,Y,Z,t)|<Kjj/J/^   u^axdydzj''  *  ^2iJIf^    J 


2, 


where  R  is  the  region  in  space  between  E  and  a  fixed  sphere  of  radius  R 
depending  only  on  X;Y,Z.  The  two  positive  constants  K, ^  Kg  depend  also 
only  on  X;  Yj  Z. 

Lemma  2.   Let  R  be  a  sphere  of  radius  R^  then  one  has 


(36)    f  /    u^dxdydz]      <  K^  /  /      |Vu|^  dxdydz 
i<R         /        ^  \~-^  r<R 
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Proof:  Let  u  be  a  solution  of  (l)  -  (5)-   Then  u   is  also  a  solution  of 
(l)  -  (5)  with  different  functions  f  and  g.   Thus  u  also  satisfies  (6)  but 
with  a  different  constant  VL.      By  lemmas  one  and  two. 


r 


(37)       |u|   <  ICK  /      |Vw|      sin  9  drdcpd9 


\    1/2 
'  + 


r^  w 


J 


K^ 


tt 


1/2 


R     c"(r) 


—       sin  0  drdcpdQ 


Applying   (6)    to   (37)    for  t   such  that  R  <  e  7t   ^    one   obtains 


|u|    <  M^t 


■(1-5/2) 


Thus  if  a  <  2  one  will  have  a  decay  for  the  solution. 
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